The Origin of the Binary Pulsar J0737-3039B 
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Evolutionary scenarios suggest that the progenitor of the new binary pulsar J0737-3039B 0, 
was a He-star with M > 2.1 — 2.3 Mq [3, 3- We show that this case implies that the binary 
must have a large (> 120 km/s) center of mass velocity. However, the location, ~ 50 pc from the 
Galactic plane, suggests that the system has, at high likelihood, a significantly smaller center of mass 
velocity and a progenitor more massive than 2.1 Mq is ruled out (at 97% c.l.). A progenitor mass 
around 1.45 Mq, involving a new previously unseen gravitational collapse, is kinematically favored. 
The low massprogenitor is consistent with the recent scintillations based velocity measurement of 
66±15 km/s [5| (and which also rules out the high mass solution at 99% c.l.) and inconsistent with 
the higher earlier estimates of 141±8.5 km/s ||. Direct proper motion measurements, that should 
be available within a year or so, should better help to distinguish between the two scenarios. 

PACS numbers: 97.60.Gb,97.80.-d,97.60.Bw,97.60.-s 



The remarkable binary system J0737-3039 0,0 is com- 
posed of two pulsars denoted A and B. We show here 
that the orbital parameters of this system and its loca- 
tion close to the galactic plane pose strong limits on the 
origin of this binary system and on the progenitor's mass 
of the younger pulsar B. 

The separation, R, (i.e., the sum of the semi-major 
axes) of the pulsars today is 8.8 x 10 10 cm and the ec- 
centricity, e, is 0.087779. Both R, and e decrease with 
time due to gravitational radiation emission and the sys- 
tem will merge in approximately 85 Myr from now. The 
periods Pa.b and their time derivatives provide upper 
limits for the life times of the pulsars: s=s 210 Myr and 
ts ~ 50 Myr. Integration backwards in time yields 
the system's parameters at birth, 50 Myr ago when B 
was born: R « 10 11 cm and e w 0.11, both only slightly 
larger than the present values. The values 210 Myr ago 
were not that different: R = 1.2 x 10 11 cm and e « 0.14. 
As these parameters did not notably evolve, our analysis 
is insensitive to the exact age of the pulsar. 

Dewi and van den Heuvel and Willems and Kalogera 
P| considered a scenario in which the progenitor star lost 
most of its envelope through interaction with its compan- 
ion A. Prior to the formation of the second pulsar, tidal 
interaction between the progenitor and the neutron star 
has led to a circular orbit. Since most of the lost mass 
could not have accreted onto the companion, it was prob- 
ably lost through a common envelope phase, at which 
point the companion J0737-3039A was spun up and its 
magnetic field was suppressed by accretion. They esti- 
mate that the progenitor mass was more massive than 
2.3 M [3 or 2.1 M respectively. This limit follows 
from standard evolutionary scenarios, leading neither to 
neutron star formation nor to core collapse, from pro- 
genitors that are less massive than 2.1 — 2.3 Mq The 
formation of the second pulsar is described according to 



the current picture by a core collapse event that involves 
a supernova and mass ejection from the system. 

We consider now the influence of mass ejection during 
the formation of the second pulsar, on the orbital motion. 
The thrust of the ejected mass, Am, gives a velocity, v cm , 
to the center of mass (CM) of the remaining system. In 
addition, the mass loss would lead to an elliptic orbit or 
even to the disruption of the system. For a spherically 
symmetric mass loss, v cm will be: 
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where vk = \jG(jriA + maj/K is the Keplerian veloc- 
ity of the two stars relative to each other, just after the 
explosion, with R being the distance between the two 
stars at that time. Within the context of J0737-3039, 
m A = 1.377(5)M Q , m B = 1.250(5)M Q , Am is, of course, 
unknown and mg, e tub + Am is the initial mass of B, 
while vk ~ 600 km/s. With ttla ~ tub, v cm would be of 
the order of Vk/2 unless Am C mg. The birth of the 
system with a low eccentricity and Am « mj requires 
for B to have had a natal kick ||| . This kick will increase 
v cm further and the above value can be considered as a 
lower limit. In other words, the nearly circular orbit to- 
day implies either a large CM velocity, roughly as given 
by Eq. E] or a small ejected mass Am <C m^. 

Fig. 1 depicts change in the CM velocity (relative to 
the unknown CM velocity prior to the explosion) for a 
system with two masses moving on a circular orbit, as- 
suming the new binary system attains the initial eccen- 
tricity of J0737-3039 (0.088 < e < 0.14). In the following 
we assume that this initial CM velocity was small and we 
approximate v cm ps Av cm . We will return to this point 
in the conclusions. Fig. 1 shows that the minimal CM 
velocity for a 2.1 Mq progenitor is 120 km/s. 

Recently Ransom et al. have estimated v cm x, the 
CM velocity of the binary on the plane of the sky, using 
the observed scintillations of the system. They find a 
rather large value: v cm ± — 141±8.5km/s (with 96.0±3.7 
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FIG. 1: The maximal kinematically allowed mass loss 
Am m « and progenitor mass msi for a given change in the 
CM velocity Av cm . The upper mass depends on the actual 
eccentricity (0.088 < e < 0.14), hence the finite width of the 
line. The hatched area is kinematically excluded. The region 
allowed by standard evolutionary scenario through the forma- 
tion of a He-star, which requires a minimum mass of 2.1 Mq 
(marked W&K Q) to 2.3 M (marked D&vdH Q) to form 
a NS through a core collapse SN, is marked as shaded trian- 
gles. The no-kick solution that results with the above range 
of eccentricities is marked by the short heavy line. 

km/s along the orbit and 103.1 ±7.7 km/s perpendicular 
to it). This value excludes the region in Fig. 1 left of a 
vertical line of ~141 km/s. However, these findings were 
questioned recently by Cole et al. [f| who suggested that 
the scintillation pattern is anisotropic. When including 
anisotropy, they find a much lower value: v cm ± = 66 ± 
15 km/s. The region in Fig. 1 to the right of the vertical 
line of ~66 km/s is consistent with this observation. 

However, the system can be additionally constrained. 
The observed distance of the system from the Galactic 
plane, z b s « 50 pc, enables us to place a statistical up- 
per limit on v cm . Stars move in a periodic motion in the 
vertical direction. For small vertical oscillations, the po- 
tential of the Galaxy is harmonic: $ = 2irGpoz 2 , where 
pa w 0.25 M Q /pc 3 is the mass density in the disk [To| . 
This gives a vertical orbital period, P z w 50 Myr. The 
typical velocity for an object at z b s is v z « 27rz D b s /P 2 . 
z bs ~ 50 pc implies then that the expectation value of 
the vertical velocity is of the order of 6 km/s. 

To quantify the probability for having a particular CM 
velocity given the observed z b s we perform Monte Carlo 
simulations that follow the formation of the system. We 
assume that star B had a given mass niBi, and that a ran- 
domly oriented kick \~kick was given to it. We also assume 
that the progenitor distribution has an initial Gaussian 
distribution in the amplitude of the vertical oscillation, 
with a width a z = 50 pc (other <7 Z <100 pc gave very 
similar results) . At the moment of formation, we assume 
it had a random phase within its vertical motion. We 
calculate the CM kick velocity v cm , and assign it a ran- 
dom direction, then integrate the vertical motion of the 
pulsar for 50 Myr using a realistic galactic potential ^3 • 

Fig. 2a depicts the probability that a system with a 
given vriBi and v^ick could find itself with 0.087 < e < 
0.14 after 50 Myr, within 50 pc of the galactic plane and 
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FIG. 2: Panel A: The probability that a binary system will 
end up within 50 pc from the galactic plane, with 0.088 < e < 
0.14 and with a transverse velocity of 66±15 km/s, given that 
50 Myrs before, the progenitor system had a circular orbit, 
that star B had a progenitor mass rriB,i, and it obtained a 
random kick velocity of size «kick- The probabilities are nor- 
malized to the most likely Mbi, ""kick- The He-star solution 
requires m,Bi larger than about 2.1M0. With a fine tuned 
i/kick, this type of a solution is ruled out at the 99% c.l. So- 
lutions with either m,Bi ~ 1-55 ± 0.2 Mq and Vkick^30 km/s, 
or, mBi ~ 1.45 ± 0.2 Mq and Wkick ~ 50 — 80 km/s, are kine- 
matically more favorable. In panel B, the constraint on the 
transverse velocity is alleviated. Here, He-star scenarios can 
be ruled out only at the 97% c.l. The top left side is kinemat- 
ically forbidden (no solution without a large «kick exists to 
counter the ejected mass), while the r.h.s. is statistically dis- 
favored since a large fJkick implies a low probability for finding 
the system near the galactic plane. Unlike Panel A, the ec- 
centricity in panel C is constrained to be e < 0.14 — showing 
that double ridged solution arises from the constraint on ec- 
centricity. Panel D is the same as panel A, except that here 
the transverse velocity is 141 ±8.5 km/s. Note that as marked 
on Panel B while the region on the upper left side - high mass 
and low kick velocity is kinematically forbidden, the region on 
the right hand side and in particular the top right hand side 
- high mass and large kick velocity is allowed kinematically it 
is just disfavored statistically. 



with a transverse velocity of 66± 15 km/s as measured by 
Coles et al. [5j. For comparison, Fig. 2b depicts similar 
simulations but without the assumption on the transverse 
velocity. Here the initial conditions are less constrained. 
Fig. 4c repeats the first panel, with the modified condi- 
tion that the e < 0.14 instead of 0.087 < e < 0.14. One 
can see that qualitatively the results do not change. Fig 
4d is a repeat of the first panel with the transverse veloc- 
ity assumed to be 141 ±8. 5 km/s, as measured by Ransom 
et al. 0, with the scintillation anisotropy neglected. 
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Without the additional information on the CM veloc- 
ity, we find that that even if Vkick is fine tuned to be 
near either 130 or 315 km/s, a system with a mass of 
iTiBi = 2.1 Mq could result with the observed configura- 
tion in only about 3% of the random realizations (as com- 
pared with the most favorable conditions having lower 
mBi and Ukick)- Other ^kick's, or higher mass systems 
are kinematically even less likely. On the other hand, 
low ejected mass solutions are favored. If we add the 
constraint that the transverse velocity is 66 ± 15 km/s, 
then a fine tuned rriBi = 2.1 Mq model can be ruled out 
at even 99% c.l. The results, are though, qualitatively 
different if we take the CM velocity on the plane of the 
sky as 141 km/s. Much larger kicks, of a few hundred 
km/s are essential now. While the "cannonical" 2.1Mq 
or higher solution is not really favorable yet, the pre- 
viously most favorable low mass (1.4M ) progenitor is 
ruled out. Instead an "intermediate" mass of ~ 1.7M 
is the most likely one. 

If we assume that the transverse measurement of 
Coles et al. which allows for the scintillation to be 
anisotropic, then at better than 99% confidence, the pro- 
genitor of star B was less massive than 2.1 Mm. This is 
just the lower limit of Willems and Kalogera Q and it is 
slightly below the lower limit of Dewi and van den Heuvel 
3] for the progenitor mass in a standard core collapse 
scenario. An inspection of Fig. 2a reveals that there is 
a kinematically favorable solution with a small mass loss 
and natal kicks ranging from to ~ 100 km/s. Since a 
small mass loss necessarily implies a small natal kick, the 
solutions with Vfci C fc>100 km/s are physically unlikely. If 
we therefore limit ourselves to Vkick^O km/s, then with- 
out fine tuning, a large fraction of the progenitor systems 
will result with the observed configuration with a progen- 
itor mass of 1.45 M Q <mBi<1.65 Mq. Thus, a mass loss 
of about 0.3 ± 0.1 Mq is most probable. The results 
are qualitatively the same even if we do not enforce that 
the resulting system has a CM velocity of 66 km/s, or if 
we replace the condition on the eccentricity to e < 0.14 
rather than 0.088 < e < 0.14. 

At a distance of 600 pc a velocity on the plane 
of the sky, v cm ±, implies a proper motion of 
0.036(v cm _L/100 km/s) "/yr. This should be compared 
with the current errors of ~ 0.04" in the position of 
the system. This comparison suggests that within a 
year or so we could obtain a direct limit on the pecu- 
liar motion. A comparison of this peculiar motion with 
the motion measured using scintillations would confirm 
the anisotropy estimates that arise from the scintillations 
measurements. The break down to the different compo- 
nents should allow us to measure the orientation of the 
orbital plane in the sky. Thus within a year or two we 
should be able to nail down the two components of the 
CM motion on the plane of the sky as well as the orien- 
tation of the orbital plane! The knowledge of the vertical 
component (relative to the galactic plane) would tell us 
whether the position of the system in the galactic plane is 
natural (if the vertical component is small) or an unlikely 



coincidence (if this velocity is large). 

We are left with two physically distinguished scenarios 
for the formation of pulsar B. In the first, the progenitor 
is a kinematically "unlikely" but theoretically plausible 
2.1 — 2.3 Mq He-star progenitor — around the minimal 
masses estimated from stellar evolution scenario 0, H|- 
In the second scenario, the progenitor is a kinematically 
favorable ~ 1.5 M young stellar core. The more proba- 
ble low mass solution requires a new type of stellar col- 
lapse. Intermediate solutions are neither statistically fa- 
vored nor do they fit any plausible theoretical scenario. 

Before turning to the implications of these two solu- 
tions we consider, first, three assumptions made in our 
analysis, (i) We have assumed that prior to the for- 
mation of the second pulsar the system was in a circu- 
lar motion. This follows from all evolutionary scenarios 
that lead to a neutron star and a small mass progenitor 
(even a 2.3M Q is a very small mass progenitor), (ii) We 
have assumed that the second mass loss was instanta- 
neous, namely, shorter than a fraction of an orbital pe- 
riod. Given that the orbital motion is of several hundred 
km/s while typical mass ejection velocities in SNe are 
higher than 10,000 km/s, this assumption is reasonable 
(for all conventional neutron star formation scenarios), 
(iii) We have assumed that the CM velocity, prior to the 
formation of the second pulsar was small. One would ex- 
pect that the system would have acquired a CM velocity 
during the formation of the first pulsar. About half of the 
system's mass was lost during this event and this should 
have resulted in some CM velocity. However, this veloc- 
ity would have been of the order of half the Keplerian 
velocity of the system at that time (see Eq. QJ. Given 
the fact that the orbital separation was much larger we 
could reasonably expect, but not prove, that this velocity 
was of order of several tens of km/s. Moreover, if the sys- 
tem would have acquired a large CM velocity in the first 
SN, there would have been an even smaller probability 
to find it in the galactic plane today. 

We turn now to the most likely scenario, the very low 
mass scenario. We imagine the same evolutionary sce- 
nario in which some time before 50 Myr, system A and 
progenitor B were in a common envelope phase and B 
lost most of its mass keeping practically just its core of 
~ 1.45 Mq. This progenitor leads to a small CM mo- 
tion and does not require any kick velocity. This solu- 
tion is kinematical preferred. However, it requires a new 
mechanism for the formation of the pulsar as He-stars of 
1.45 Mq do not collapse to form neutron stars 0, Il3| . 
The observation that the progenitor mass is very close 
to the Chandrasekhar mass leads us to conjecture that 
the process involves the collapse of a supercritical white 
dwarf. For example, the progenitor may have been a de- 
generate bare core just above the critical Chandrasekhar 
mass, supported by the extra thermal pressure against 
collapse. As it cooled, the additional support was lost 
and the core collapsed to form a neutron star. A second 
possibility is that it was formed just below the Chan- 
drasekhar mass, and as it cooled, neutronization at the 
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core increased the baryon to electron ratio, and with it 
reduced the Chandrasekhar mass until the progenitor be- 
came unstable and collapsed. Note that the object must 
have been composed of O-Ne-Mg as a collapsing CO core 
would have carbon-detonated and it would have exploded 
completely forming a type I SNe and leaving no remnant. 
This solution clearly requires a new type of formation 
scenario for neutron stars. 

This new solution passes an immediate non trivial 
test. With a small mass loss, only small kick veloc- 
ities are possible and in fact we can estimate in this 
case (see Fig. 2) the mass loss needed to obtain the 
initial eccentricity e w 0.11. We find (while conserva- 
tively assuming that the collapse took place 50+5°° Myr 
ago) Am = e(l + q)m B = 0.28 ± 0.07 M and a pro- 
genitor mass of 1.53 ± 0.07 M , which is, indeed, just 
above the Chandrasekhar limit. Intriguingly, some mass 
loss, in the form of v losses of a few times 10 53 ergs, 
must take place. The estimated Am corresponds to 
E v w Amc 2 m 4.2 x 10 53 ergs, which is in the right range. 
Of course, if some mass is ejected as well, E v will be 
smaller and a small kick could arise, but the total mass- 
energy lost will be the same. The consistency of this 
mass and energy loss with the previous physical picture 
increases our belief in this new and unusual scenario. 

The other, higher mass scenario of 2.1 — 2.3M Q (which 
is still on the lowest end of the evolutions scenarios) is 
the most likely from a stellar evolution point of view. 
However, as we have seen it is quite unlikely statistically. 
Specifically, it implies a post-kick CM velocity of at least 



125 to 150 km/s, which was within a few degrees from the 
Galactic plane for it to find itself in the plane today. To 
be consistent with the low v cm ± , it should also be point- 
ing within about 25° from us, or in the opposite direction. 
In other words, it should have a minimum velocity of at 
least 100 km/s in the direction of the line of sight, if pro- 
genitor B was a 2.1M star. More massive progenitors 
imply even larger velocities. Of course it cannot be ruled 
out. We might be observing a very improbable system. 
However, if this is the case, the "unlikeliness" of the sys- 
tem should be taken into consideration when performing 
population synthesis and estimating neutron star merg- 
ers rate based on this observation. Specifically, with a 
very large CM velocity, as such a system must have, it 
is very unlikely to find it in the Galactic plane. Most 
such systems should exist high above the Galactic disk. 
Anticipating this, Narayan, Piran & Shemi 0] consid- 
ered this effect by taking an effective scale hight of 5kpc 
for the Galactic binary NS systems. This value was used 
later by others. It may be that even this high value is 
too low. This in turn might boost up the estimated rates 
of NS merger by a significant factor. 
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